arXiv:1505.05728vl [gr-qc] 21 May 2015 


Isofrequency pairing of spinning particles in Schwarzschild-de Sitter spacetime 
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It has been established in Schwarzschild spacetime (and more generally in Kerr spacetime) that 
pairs of geometrically different timelike geodesics with the same radial and azimuthal frequencies 
exist in the strong field regime. The occurrence of this socalled isofrequency pairing is of relevance 
in view of gravitational-wave observations. In this paper we generalize the results on isofrequency 
pairing in two directions. Firstly, we allow for a (positive) cosmological constant, i.e., we replace 
the Schwarzschild spacetime with the Schwarzschild-de Sitter spacetime. Secondly, we consider not 
only spinless test-particles (i.e., timelike geodesics) but also test-particles with spin. In the latter 
case we restrict to the case that the motion is in the equatorial plane with the spin perpendicular 
to this plane. We find that the cosmological constant as well as the spin have distinct impacts on 
the description of bound motion in the frequency domain. 
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1. INTRODUCTION 

The occurrence of isofrequency pairs in the strong field 
regime of the Schwarschild spacetime was noticed only re¬ 
cently by Barack and Sago [J. It says that bound orbits, 
for which both the radial and the azimuthal motion are 
periodic, are not uniquely characterized by their frequen¬ 
cies. More precisely, there exist geometrically distinct 
orbits in the strong field regime that possess the same 
frequency pairs. 

At first, this degeneracy feature may not be much of 
a surprise. After all, it is known from Newtonian Me¬ 
chanics that the frequencies of the Kepler ellipses are 
all degenerate, i.e. the radial and azimuthal frequencies 
have the same value. This is the reason why the orbits 
in Newtonian physics are closed. 

When general relativistic effects are considered, 
though, one major difference to Newtonian physics is 
the periastron shift of bound orbits which are no longer 
closed. This manifests itself in the non-degeneracy of 
the frequencies. In the Schwarzschild spacetime we have 
two independent orbital frequencies, for the radial and 
for the azimuthal motion. It was long thought that these 
two frequencies provide another unique parametrization 
of the orbits, as an alternative to the ones already known, 
such as the energy and the angular momentum or the pe¬ 
riastron and the apastron. However, Barack and Sago [I| 
showed that in the strong field of Schwarzschild space- 
time, i.e. in the highly relativistic regime, there ex¬ 
ist pairs of timelike geodesics which are described by 
the same frequencies. In a follow-up study, Warburton, 
Barack and Sago [H generalized the isofrequency pairing 
to the Kerr geometry. In contrast to the Schwarzschild 
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case, timelike geodesics in the Kerr spacetime have three 
degrees of freedom and in their study Warburton et al. 
found even triperiodic partners. As outlined in Refs. 
Hi, the occurrence of isofrequency pairing is of rele¬ 
vance in view of gravitational wave analysis because it 
implies that, for the case of an Extreme Mass Ratio In¬ 
spiral (EMRI), from the observation of the fundamental 
frequencies one cannot uniquely determine the shape of 
the orbit. Shortly after Refs. HiUiad appeared, Shay- 
matov, Atamurotov and Ahmedov [Sj investigated the in¬ 
fluence of a magnetic field on the orbital frequencies of a 
charged particle moving in Schwarzschild spacetime and 
found that the region where isofrequency pairing occurs 
shrinks for high values of the magnetic field. 

In this paper we generalize the results on Schwarzschild 
isofrequency pairing in two different directions. Eirstly, 
we consider the Schwarzschild-de Sitter spacetime, which 
is the unique spherically symmetric and static solution to 
Einstein’s vacuum field equations with a positive cosmo¬ 
logical constant. Secondly, we consider not only spin¬ 
less but also spinning test paricles. We restrict to the 
case that the particle moves in the equatorial plane, 
with the spin perpendicular to this plane. Properties 
of the Schwarzschild-de Sitter spacetime have been dis¬ 
cussed, e.g., in d-Q, and the motion of non-spinning test- 
particles in this spacetime has been investigated, e.g., 
in HI- Eor spinning particles we have to consider the 
equations of Mathisson-Papapetrou-Dixon (Tol - [l^ . As 
this set of equations is not closed in the sense that there 
are less equations than variables, an additional condition 
is needed, usually referred to as a spin supplementary 
condition (SSC). Choosing an SSC is associated with fix¬ 
ing a reference frame the motion is observed in. The 
choice of the SSC depends on the question one wants to 
investigate [l3l-l^. We use the Tulzcyjew condition (T 
SSC) [2ll| which corresponds to the zero 3-momentum 
frame. The dynamics of spinning test-particles has been 
worked out in Lagrangian and Hamiltonian formalisms 
of which the latter appears to be more complicated Hi- 
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[^ . For the motion of spinning test-particles in the 
Schwarzschild or Kerr spacetime we refer, e.g., to [25l - 
and for the case that also the quadrupole moment 
of the test-particle is taken into account to [1^. Several 
papers have been devoted to the question of whether the 
motion is chaotic, which also has a great impact on the 
analysis of possible gravitational wave signals [29l - l3^ . 

The paper is organized as follows. In the second section 
we recall the equations of motion for non-spinning and 
spinning test-particles in Schwarzschild-de Sitter space- 
time. In particular we specify the region of bound mo¬ 
tion. After that we analyze the frequencies of bound 
motion and find the domain where isofrequent orbits ex¬ 
ist in section three. Lastly, we summarize our results and 
give an outlook for future work and applications. 


2. MOTION IN SCHWARZSCHILD DE SITTER 
SPACETIME 

2.1. Non-spinning particles 

Schwarzschild-de Sitter spacetime is the unique spheri¬ 
cally symmetric vacuum solution to Einstein’s field equa¬ 
tions including a positive cosmological constant, A > 0, 
see e.g. Rindler [^. For a negative cosmological con¬ 
stant it is called Schwarzschild-anti-de Sitter spacetime, 
but here we are interested only in the case A > 0. The 
metric reads 


ds^ = —/ (r) di^ + / (r) ^ dr^ (dO^ + sin {9)^ dtp'^'j (1) 
in spherical coordinates with 


r / ^ 1 2M a 2 

/ r = 1--rA 

r S 


( 2 ) 


Here M corresponds to the mass of the gravitating body 
and the units are chosen such that c = 1 and G = 1. We 
consider only the region r > 0. For 0 < A < (3M)“^ , 
the function /(r) has two zeros, indicating two horizons, 
at radii thi and rH 2 with 2M < thi < 3M < rH 2 ', 
the metric is static in the region rni < r < rH 2 where 
/(r) > 0. For A > (3M)“^ there is no horizon and no 
static region. 

The dynamics of massive non-spinning test-particles is 
determined by the geodesic equation for timelike curves. 
Since the spacetime is spherically symmetric the mo¬ 
tion takes place in a plane, which we choose to be the 
equatorial plane 6 = ^. Then the equations of motion, 
parametrized by proper time r, for a test-particle of mass 
m read 


y / dt\'^ 


(3) 


f{r)^ ’ 

a / dr \ 

, ^J^ 



= H^-f{r) - 

(4) 

2(d(j>\^ 

= ^1 
r4 ’ 

(5) 


with H and Jz being two constants of motion to be in¬ 
terpreted as the energy and the angular momentum, re¬ 
spectively. The particle’s four-momentum is 
and satisfies the mass-shell condition rn? = —p^p^. From 
hereon we rescale H Hm and Jz >—t Jz'm which is tan¬ 
tamount to setting m = I in ©, m and ([5]). Then H is 
dimensionless while Jz has the dimension of a length. 

In the following we are interested only in bound mo¬ 
tion, i.e., in orbits that have two turning points, an apas- 
tron Ta and a periastron Tp, where dr/dr = 0. From (|3]) 
it is obvious that turning points can exist only at radial 
coordinates where /(r) > 0. This implies that bound or¬ 
bits are confined to the region between the two horizons. 
In particular they do not exist if A > (3M)~^. 

We rewrite dU as 




A 


( 6 ) 


where 


PsW =r 5 - ((l-ij 2 )__ j; 


6M 


- - J^r 


6M 


Ji 


(7) 


As we assume A > 0, the region where P 5 {r) < 0 is for¬ 
bidden by (|ni). The number of zeros of P^ir) determines 
the types of motion possible for the corresponding values 
of H and Jz- For positive A there can be at most four 
positive real zeros as can be derived with the rule of signs 
by Descartes. Bound orbits are allowed only when the 
polynomial has precisely four positive zeros. For other 
values of PI and Jz there will be escape or terminating 
orbits or no motion at all. 

From now on we rescale r i—> rM, Jz i—> JzM and 
A I—>■ AM~‘^ so that these quantities are dimensionless. 
(Recall that PI already was dimensionless.) This is tan¬ 
tamount to setting M = 1 in ([7]). Since we must have 
four positive real zeros in order to consider bound motion 
we can rewrite equation o as 


P^(r) = {r - ro)(r - ri)(r - rp)(r - ra){r - r 2 ) (8) 

with To < 0 < ri < Tp < Ta < r 2 . Bound motion 
exists only between Vp and r^. Also, the motion is fully 
described by the parameters ra and rp, because P[ and 
Jz can be expressed in terms of ra and rp, 

(t^g -ra + 3) (f rg - rp -F 3) (r„ rp) 

rl{rp - 2) + r'^ {ra - 2) - 2rarp 
rlr^jS - ^rgrpjra + rp)) 

- 2) -b r-2 {ra - 2) - 2rarp ' 

Therewith we can find all the zeros of © expressed in 
terms of rp and Comparing the coefficients of the 
two equations © and © yields expressions for rg and 
r 2 dependent on {ri,ra,rp) and the equation 

P3{ri) := Arlr^rf + Arlrl{ra + rp)rl 

+Jl{&rp-2>ra{rp-2))ri+&Jlrarp = 0 (II) 
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so that we only have to solve a cubic equation for ri. 

In the {H, J^)—plane, the boundaries of bound motion 
are determinded by the merging of two zeros or by the 
condition that the orbit becomes unbound. Physically, a 
merging corresponds to a stable circular orbit. 

A merging ri = rp ^ 7 ^ r 2 corresponds to an unstable 

circular orbit at rp and a homoclinic orbit from rp to 
Ta and back to Cp, while a merging ri ^ Vp ^ Va = ^2 
corresponds to an unstable circular orbit at ra and a 
homoclinic orbit from Va to Vp and back to Tq. In the 
case that ri = Vp ^ Va = r 2 we have unstable circular 
orbits at rp and at and heteroclinic orbits from to 
Tp and from rp to r^. 

Another parametrization of bound orbits is given by 
the semi-latus rectum p and the eccentricity e. The rela¬ 
tions to the periastron rp and the apastron ra are 


ra = 



P 

1 -I- e ' 


( 12 ) 


We choose this parametrization since then the analysis 
of the properties of a test particle’s motion is more con¬ 
venient. 


e 



Figure 1: The figure presents the regions of bound motion in 
the (p, e) - plane for different values of A. The dashed line is 
the separatrix in the Schwarzschild case, A = 0, given by the 
equation p — 6 -|- 2e. The remaining lines are the separatrices 
for a cosmological constant of A = 0.0001 (upper dotted line), 
A = 0.0005 (solid line), and A = 0.0006 (lower dotted line). 


In Fig. [T] the region of bound motion is shown in a 
(p, e)-diagram for both Schwarzschild and Schwarzschild- 
de Sitter spacetimes. We first consider the Schwarzschild 
case for which the general features of the domain of 
bound orbits have been discussed already by Cutler, Ken- 
nefick and Poisson [sj, see also [ij and Hi- In this 
case the domain of bound orbits, corresponding to the re¬ 
gion that lies to the right of the dashed line, is infinitely 
large. This dashed line, which is given by the equation 


p = 6 + 2 e, is often called the separatrix because it sep¬ 
arates the region of bound orbits from the region of un¬ 
bound orbits. On the separatrix, each point corresponds 
to a homoclinic orbit from a radius ri = rp to a radius 
ra > rp and back to ri = rp. Such a homoclinic orbit 
has the same constants of motion H and Jz as the un¬ 
stable circular orbit at ri = rp which it asymptotically 
approaches. Bound orbits near the separatrix are period¬ 
ically “zooming out” to the apastron with many “whirls” 
near the periastron in between, see Levin, O’Reilly and 
Copeland 
cussed in 


Ml such zoom-whirl orbits have also been dis- 
42-1^. The upper boundary curve of the re¬ 


gion of bound orbits corresponds to unbound (parabolic- 
type) orbits with e = 1 , while the lower boundary cor¬ 
responds to stable circular orbits with e = 0. The lower 
left-hand corner of the region of bound orbits corresponds 
to the innermost stable circular orbit (ISCO) at r = 6 . 


By contrast, in the Schwarzschild-de Sitter spacetime 
the region of bound orbits in the (p, e)—plane is finite. 
The shape is triangle-like with its tip at an eccentricity 
Cmax < 1- This demonstrates that, in this picture, the 
transition from bound orbits to unbound orbits (e > 1 ) 
is not continuous. In analogy to the Schwarzschild case, 
the two sides of the triangle are called the separatrices. 
They correspond to homoclinic orbits with ri = rp and 
ra = ^ 2 , respectively. Each homoclinic orbit has the same 
values for energy and angular momentum as the unsta¬ 
ble circular orbit that it approaches asymptotically. We 
have already mentioned that in the Schwarzschild space- 
time bound orbits near the separatrix show a zoom-whirl 
behavior, with the “whirling” taking place in the strong- 
held regime. In the Schwarzschild-de Sitter spacetime 
with a small positive cosmological constant, the second 
separatrix gives rise to zoom-whirl orbits that “whirl” 
near an apastron far away from the center and period¬ 
ically “zoom in” to a periastron. The two separatrices 
intersect at the tip of the triangle where we have si¬ 
multaneously ri = rp and = r 2 . This gives rise 
to a heteroclinic orbit, i.e., to an orbit that connects 
two different unstable circular orbits. Such an orbit can 
be described by a unique pair of values for (p, e). The 
lower boundary curve of the region of bound orbits cor¬ 
responds to stable circular orbits. The two intersection 
points of the separatrices with the horizontal axis corre¬ 
spond to the innermost stable circular orbit (ISCO) and 
the outermost stable circular orbit (OSCO). While the 
ISCO is also present in the Schwarzschild spacetime and 
merely gets shifted away from the center with increasing 
A, the OSCO is only existent in the Schwarzschild-de Sit¬ 
ter spacetime. For A —>■ 0 the OSCO approaches infinity. 
For A —>• Acrit = 4/5625 the ISCO and the OSCO merge 
into one circular orbit, cf. Q. Put into mathematical 
language, this happens if all four positive zeros of © co¬ 
incide. For A > Acrit bound orbits do not exist. Let us 
remark here that Acrit is much larger than the physically 
expected value of the cosmological constant. Observa¬ 
tions show evidence for a A Ri [s^. Even for 

a supermassive black hole with M r; lO^^km this corre- 
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spends, in our geometrized units, to A r; 10 i.e., to a appropriate adaptation, 

value that is much smaller than Acrit- 


2.2. Spinning particles in the equatorial plane 


Now we turn on the spin and investigate the motion of 
a spinning test-particle which is no longer geodesic. The 
corresponding equations of motion are the well-known 
Mathisson-Papapetrou-Dixon equations [iMl, 


_ 1 pAi okX 

P — ’ 

(13) 

= pu’' - p’'u>^. 

(14) 


Here is the curvature tensor of the underlying 

spacetime, the dot denotes differentiation with respect 
to proper time r along the worldline of the particle, 
is the four-momentum and is the four-velocity of the 
particle, and is the antisymmetric spin tensor. As 
already mentioned in the introduction, the set (fO]) - (n3 
is underdetermined; a so-called supplementary spin con¬ 
dition (SSC) must be chosen in order to get a well-posed 
initial-value problem. In this work the Tulzcyjew condi¬ 
tion [2l| is chosen. 


= 0 . 


(15) 


Then the relation between p^ and u'^ can be given explic¬ 
itly: Using the conserved quantity = —PfiP^, intro¬ 


ducing the normalized momentum vector P^ = ^ and 
renormalizing the four-velocity such that = —1 al¬ 
lows to derive the relation 




2 [m? + • 


(16) 


Therewith, it is now possible to derive the equations of 
motion for a spinning test-particle in Schwarzschild-de 
Sitter spacetime. However, here we will restrict ourselves 
to the special case of a particle moving in the equatorial 
plane, with the spin vector perpendicular to this plane. 
Then we can characterize the spin by the scalar constant 
of motion S', defined by 

S^ = S^S^ , S'^ = (17) 


with the totally antisymmetric Levi-Civita tensor field, 
and the property that S is positive if the spin is paral¬ 
lel to the orbital angular momentum and negative if it is 
anti-parallel. As in the spinless case, we have a conserved 
energy H and a conserved angular momentum which 
we rescale according to H Hm and J^^m. In 

addition, we now also rescale the spin, S i—> sm. The re¬ 
sulting equations of motion can be taken from , after 


^ H+^sJ, 

Kdr) Ils{r)'Es{r)f{r) ’ 

Kdr) ns(r)Es(r-) ’ 

_ (■/. - Hs) (l - -4^5^) 

Vdr/ ns(r)Es(r)2r2 


( 18 ) 

(19) 

( 20 ) 


Here /(r) is defined by ([2|), the prime denotes derivative 
with respect to r, and 


n^(r) 

^s{r) 

7?s(r) 



I - 



(J, - Hsf 
2r^Es(r')3 


( 21 ) 

( 22 ) 



fir) ^ 

S.(r)2 + 


(Jz - Hsf ' 


(23) 


Again we are interested in bound motion, so we require 
two turning points and Vp where dr/dr = 0 which 
corresponds to i?s(r) = 0. We rescale, as before, r >->• 
rM, Jz I—>■ JzM, A I—>■ AM~^ and now also s i—>■ sM. 
Note that for test particle motion our dimensionless spin 
parameter s necessarily satisfies the condition — 1 < s < 
1, see [ 13 . Moreover, we substitute 


L = Jz-Hs 


(24) 


for mathematical convenience. Then Rs (r) can be rewrit¬ 
ten as 


Rs{r) 


3r^ 


Pfr) 


(25) 


where 


Pgir) = r® -|- ar^ -|- 5r® -|- cr® -1- dr'^ er^ gr h (26) 
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with 


a 

b 


c 


d 


e 


9 

h 


(i/(1 + l-s^) + Ais) +1-1/2 g 

t (1 + ^^)' ^ 


+2 

2(^{L + Hsf - 2L(2L + Hs) 

f ^ f (l + fs^)' 

4s^(1 + a^ 
s2 ^(L + ifs)^ — 

2s4 

t (1 + ^^)" 


(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


The number of zeros of Pgir) determines the types of 
motion possible in the corresponding spacetime. Finding 
the maximum number of positive real zeros in this case is 
not so easy, though, since the signs of a and d are unclear. 
(The sign of g does not matter because the signs of e and 
h are already different.) In any case, from Descartes’ rule 
of signs we hnd that the number of positive real zeros can 
only be 0, 2, 4 or 6. Again, in order to consider bound 
motion we must have at least four positive real zeros. We 
argue that in this case (IMll can be written as 


Pg (r) = (r - r^){r - rp){r - r„)(r - ri) 
x(r - rp){r - ra){r - r 2 ){r - r/)(r - fj) 

with Tg, < < Tq, < 0 < ri < Tp < Tq < r 2 and 

r/, fj being some non-real complex zero and its conju¬ 
gate, respectively. The knowledge of a non-real zero was 
obtained by checking the signs of the coefficients in the 
relevant parameter space. To exclude the case with six 
positive real zeros, we observe that in the relevant range 
of parameters, 0 < A < 1/9 and — 1 < s < 1, only three 
mergers of zeros exist. Therefore we conclude to have 
one pair of complex zeroes. 

Bound motion exists only between Vp and r^. Since 
also for the spinning particle the motion is fully de¬ 
scribed by the parameters and rp, we find expressions 
for H{ra,rp) and L{ra,rp) by setting Rs{ra) = 0 and 
Rs{rp) = 0. With the help of (IT^ the expressions are 
converted to H{p,e) and L{p,e). The next task is to 
find the boundaries of the region of bound motion in the 
(p, e)-plane, i.e. the separatrices. Again, the boundaries 
are given by the merging of two zeros, ri = rp or/and 
T'p = ’’a or/and Va = r^. Unfortunately, the explicit ex¬ 
pressions for the zeros are not as easily found as in the 
non-spinning case. However, there is a way around: We 


simply exploit the fact that the values of H and L for the 
homoclinic orbits are identical to the ones for the unsta¬ 
ble circular orbits which they approach asymptotically. 
Having access to the values of H and L for circular orbits 
by solving Rs{r) = 0 and i?((r) = 0, the only computa¬ 
tion needed is to find the intersection point of the lines 
of constant H and constant L in the (p, e)-plane. In this 
way we obtain the values for p and e of the homoclinic 
orbits corresponding to the separatrices. 



Figure 2: The figure presents the regions of bound motion in 
the (p, e)-plane for different values of s for A = 0 as well as 
A = 0.0005. The straight lines on the left correspond to A = 0 
and the triangles on the right correspond to A = 0.0005. In 
either case, the (green) dashed line corresponds to s = 0.1, 
the black solid line to s = 0, and the (blue) dotted line to 
s = -0.1. 

In Fig. [2] the region of bound motion is shown in a 
(p, e)-diagram for both a spinning and a non-spinning 
particle moving in Schwarzschild and Schwarzschild-de- 
Sitter spacetime. Here, we have fixed the cosmological 
constant either to A = 0 or to A = 0.0005 and varied the 
value of the spin parameter. 

From Fig. [5] we read that for a spinning particle in 
Schwarzschild spacetime the region of bound motion is 
inhnitely large as it is for geodesic motion. The well- 
known shift of the ISCO due to the spin is visible, such 
that for positive spin it is moved inwards and for negative 
spin outwards, cf. and also [11,143. This reflects the 
coupling of the particle’s spin to its orbital angular mo¬ 
mentum. If they are parallel to each other the resulting 
force is repulsive, while it is attractive if they are an¬ 
tiparallel [391 . Notice that the upper boundary is given 
by e = 1 which corresponds to parabolic orbits and marks 
the transition from bound motion to unbound orbits. By 
and large, the general shape of the separatrices resem¬ 
bles the one for non-spinning particles in Schwarzschild 
spacetime. 

When a positive cosmological constant is considered, 
the value of the maximal eccentricity becomes smaller for 
negative and larger for positive spin. Correspondingly, 
the critical value of A is also shifted: for positive spin 
Acrit is bigger than for the spinless case and for negative 
spin it is smaller. The dependence of Acrit on the spin is 
shown in Fig. [31 
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Figure 3: The figure presents the dependence of the critical 
value of A on the spin parameter. It monotonically increases 
with the spin. The boundaries s = — 1 and s = 1 are due to 
physical restrictions. 


It reveals that for A < Acrit(s = 0) = 4/5625 bound 
motion is possible for all positive spin values but not 
for all negative spin values. This is the reason why for 
our particular choice of A = 0.0005 it is not possible 
to have bound orbits with spins that are smaller than 
Ri — 0.5. As soon as the chosen value of A drops below 
Acrit(s = —1) ~ 0.0004 bound motion is possible for all 
spin values — 1 < s < 1. 

The influence of the spin on the ISCO in the 
Schwarzschild-de Sitter spacetime is similar as in the 
Schwarzschild spacetime. Again we see that for positive 
spin parameters the ISCO gets shifted towards the center 
and for negative spin parameters away from the center. 
For the OSCO it is the other way around. By taking into 
account these two characteristics as well as the shift of 
the maximal eccentricity, it can be immediately seen that 
the region of bound orbits becomes smaller for negative 
spin and gets larger for positive spin. 

One might think that a large positive spin is somehow 
able to destroy the existence of the heteroclinic orbit sit¬ 
ting at maximal eccentricity. However, even if the spin 
is chosen to have its maximal value of 1 the shape of 
the region of bound motion does not change. The tri¬ 
angle survives and with it the heteroclinic orbit. If A 
approaches zero the maximal eccentricity goes to 1 and 
the OSCO to infinity, for any spin value. In this case 
the separatrix resembles the one of Schwarzschild, only 
shifted closer to the center for s > 0 and farther away 
from the center for s < 0. 

In the following section we change from the 
parametrization of bound orbits in terms of p and e to 
the frequency domain. This will allow us to discuss the 
phenomenon of isofrequency pairing. 


3. ISOFREQUENCY PAIRING 

As mentioned before, Barack and Sago [l| observed 
that for bound orbits of spinless test-particles in the 


Schwarzschild spacetime the transformation from the 
constants of motion H and Jz to the radial and azimuthal 
frequencies becomes degenerate in the highly relativistic 
regime. In particular, it was demonstrated that orbits 
exist with different eccentricities but with the same fre¬ 
quency pair. 

Here, we show how the picture of this degeneracy phe¬ 
nomenon changes if a cosmological constant is turned on 
and if the particle’s spin is incorporated. Generally, the 
frequencies are defined as [s^ 




27r 
7^ ’ 

j. r 



with 


Tr = 2 



dt 

dr 


dr = 2 



dt 

dr 




(35) 


(36) 

(37) 


In correspondence with our choice of units in the previous 
sections we rescale t and analogously —>■ 


3.1. Non-spinning particles 

Again, we start with the geodesic motion. From equa¬ 
tions (0)-® we obtain 



r^dr 

/(r)\/rP5(r) ’ 


dr 

yJrP^{r) ’ 


(38) 

(39) 


with /(r) from ® and P^{r) from d?]). As the polyno¬ 
mial rP^ (r) under the root in the denominator of the 
integrand is of order 6, the integral is of hyperelliptic 
type, which cannot be integrated in terms of elementary 
functions. 

If we use, as before, the (p, e) parametrization of the 
bound orbits, we may substitute the integration variable 
r in (1551) and (1551) according to 


P 


1 + e cos X 


(40) 


where the new integration variable x is the relativis¬ 
tic anomaly. Then the boundary values of the integral 
change to 0 and tt. 

Since we want to compare the frequencies of different 
orbits, we choose in analogy to Warburton et al. Q the 
(Q^, e) parametrization for our analysis of bound orbits, 
which we are allowed to because monotonically de¬ 
creases with p if the eccentricity is held fixed. In order to 
do this we have to deal with several obstacles. First, we 
cannot analytically invert the integral to obtain p as 
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a function of e and To circumvent this hindrance the 
value for e is fixed in the integral for the frequency, so 
that the value for p can be computed using a root-finding 
method for any allowed In this way we obtain the 
value for p for any given e and that is to say, we nu¬ 
merically acquire a function p(e,n^). Hence, the radial 
frequency can also be written as a function Hr(e, and 
we are able to plot contour lines for constant into an 
e)-diagram. 

Secondly, we encounter problems close to the separa- 
trices. Both and A(j> diverge at the separatrices mak¬ 
ing it numerically challenging to perform the computa¬ 
tions for the frequencies. Luckily, there exists an approx¬ 
imation scheme for hyperelliptic integrals developed by 
Sochnev based on the approximation of irrational 
numbers. We rewrite the integrals of (155)) and (1551) with 
the substitution 


{ra + rp) + {ra - rp)x 


to obtain 


Tr =A[r^ 




Vt{x)dx 
y/Vr{x) ’ 


Acj)=B{rp,ra,A) f ^== , 

J-l yVr{x) 


(41) 


(42) 

(43) 


where Vt{x) is a rational function whose denominator has 
no zeros in the integration interval and 

Vr{x) = (l-a;^)(l-|-fcia:)(l-|-/c 2 x)(l-|-fc 3 a;)(l-|-fc 4 a:) (44) 


with 


fci 

k2 

ks 

^4 


Tg - rp 
rg +rp' 

_ rg - rp _ 

(rg - To) + {Xp - To) ’ 
_ rg - Tp _ 

{rg - ri) + {Xp - Xi) ’ 

_ rg - Xp _ 

{r2 - rg) + {x2 - Xp) ’ 


which satisfy 0 < fc 2 < fci < fca < 1 and — 1 < ^4 < 0 . 
If we get close to the separatrices either or —k^ ap¬ 
proaches 1. For each of these two cases Sochnev’s method 
allows to approximate the integrals (1421) and (1451) in terms 
of elementary integrals, see the Appendix. Consequently, 
the frequencies close to the separatrix can be approxi¬ 
mately expressed in terms of elementary functions. Now 
we have the tools we need for analyzing the behavior of 
the frequencies in the region of bound motion. 

In Fig. 0] the region of bound orbits is shown in the 
(H^, e)-plane. It is the shape of this region that is strik¬ 
ing. It looks no longer like a triangle, as in the {p, e)- 
representation, but more like a trapezoid. The tip of the 
triangle is stretched out to a straight line at Cmax- In the 
(p, e)-representation the heteroclinic orbit corresponds to 


e 



Figure 4: The figure depicts the phenomenon of isofrequency 
pairing for bound orbits in the (fl<^, e)-plane for a fixed 
A = 0.0005. The thick (red) boundary lines correspond to 
the separatrices and confine the region of bound motion. The 
(blue) solid lines inside this region coorespond to constant val¬ 
ues of Qr- The (green) dashed lines are the singular curves, 
i.e. the locus where the Jacobian determinant of the transfor¬ 
mation from (p, e) to (Ht.,H,^) vanishes. Isofrequency pairing 
occurs in two domains each of which is bounded by a sepa¬ 
ratrix and a curve of Circular Orbit Duals (COD). The COD 
curves are shown by the dotted black lines. 


a single point - the tip of the triangle - since it has a 
uniquely defined pair (r^, Xg). By contrast, the azimuthal 
frequency is not uniquely defined for the heteroclinic or¬ 
bit. Note that the original definition of the frequencies 
is valid only within the region of bound orbits. On the 
boundaries, that is at the separatrices, the frequencies 
are defined only by a continuous extension which assigns 
unique frequencies to the homoclinic orbits. However, in 
the case of the heteroclinic orbit the value of depends 
on how the orbit is approached. This is the reason why 
the heteroclinic orbit is stretched out to a straight line 
in the (H,^, e)-diagram. A bound orbit that lies close to 
this line could be called a “whirl-whirl orbit” because it 
periodically changes between a large number of whirls 
near its apastron and a large number of whirls near its 
periastron. 

The main property we are interested in is the isofre¬ 
quency pairing. This phenomenon is easily seen in the 
diagram close to the two separatrices. Looking at a con¬ 
tour line for = constant near one of the separatrices, 
we see that it has two intersection points with a contour 
line for = constant, i.e., with a vertical line in this 
diagram. Since these intersection points correspond to 
orbits of different eccentricities, we conclude that there 
exist two geometrically distinct orbits with the same pair 
of frequencies. In mathematical terms this means that 
the transformation from the frequencies ( 12 ^, to (p, e) 
is not one-to-one. In order to prove this degeneracy it is 
sufficient to show that the Jacobi determinant 


9(Hr,H0) 

d {p, e) 


( 45 ) 


















becomes singular somewhere within the region of bound 
orbits. In Fig. 2] these singular points can be found as 
the points where the tangents to the contour lines of 
rir = constant become vertical. This happens along the 


two (green) dashed curves in Fig. |4] which are called the 
singular curves. To verify that J does have two zeros 
close to e = 0, one may perform a Taylor expansion of J 
about e = 0 up to first order, 


J = 


ePio {p) 


- l) {jP^ - P + 2) (f (4p - 15)p3 - p + 9) 


3/2 


+ 0{e^) 


with 

Pio (p) = _ 50 aV + 30AV - SlSA^p^ + 9A(80A - l)p® - dSAp^ + 918Ap'‘ - 224lAp^ + I08p^ - 1053p + 2322. 


Numerically one finds that the tenth order polynomial 
Pio ip) has precisely two positive real zeros lying withing 
the allowed range of p values for bound motion, for all 
0 < A < 1/9. 

The isofrequency pairs lie on opposite sides of one of 
the singular curves, i.e., each orbit that is located be¬ 
tween a separatrix and a singular curve has a partner or¬ 
bit on the other side of the corresponding singular curve 
which is geometrically different but has the same pair of 
frequencies. The regions where isofrequency pairing oc¬ 
curs are bounded by the socalled “Circular Orbit Dual” 
(COD) curves which are represented by the black dashed 
lines. A point on a COD curve corresponds to an orbit 
that has the same frequencies as a stable circular orbit 
situated between one of the singular curves and the cor¬ 
responding separatrix. 

In comparison to the Schwarzschild spacetime, the 
most interesting new feature of isofrequency pairing in 
the Schwarzschild-de Sitter spacetime is in the fact that 
this phenomenon now occurs not only in the highly rela¬ 
tivistic regime close to the center but also in a region far 
away from the center. 

While we could read from Fig. [T]how the ISCO radius, 
the OSCO radius and the maximal eccentricity depend on 
A, Fig.|4]gives us information on the frequencies of bound 
motion. A quite interesting property is the existence of a 
maximal radial frequency. From Fig. |4]we read that the 
lines of constant have the shape of semicircles which 
become smaller in the center of the diagram. The value 
of Dr varies between = 0 on the separatrices and a 
maximal value when the semicircle is contracted to just 
one single point. Fig. [S] shows a parametric plot of the 
maximal value of fir against p, where the parameter is the 
cosmological constant. As fir takes its maximal value on 
the axis e = 0, i.e., for the limiting case of a circular orbit, 
p is the same as the radius coordinate. We see that the 
position where fir takes its maximal value varies between 
p = 8 for A = 0 and p = 7.5 for A = Acrit = 4/5625. 

Moreover, we can read from Fig. [5] that the presence of 
a positive cosmological constant introduces another kind 
of degeneracy. Notice that fl™^^ decreases slightly when 





Figure 5: The figure shows how the value of and the 

position where this value is taken depend on A. The (red) 
dot marks the Schwarzschild case with A = 0. Along the 
solid black line A increases up to its maximal value of 1/9 
where the bound orbits vanish. 


A is slightly increased. Transferring this change to Fig.|4] 
we can think of the lines of constant fir getting shifted to¬ 
wards the center. Imagine that we have a frequency pair 
close to the right separatrix. Choose the outer contour 
line of fir = constant and let some vertical line of con¬ 
stant fl^ intersect it twice and mark these two points. As 
described above, they are two orbits with the same fre¬ 
quencies. If the cosmological constant is changed a little 
bit the line of constant fir gets shifted either to the right 
if A is reduced or to the left if A is amplified. However, 
we would still have two intersection points with the ver¬ 
tical line representing the azimuthal frequency we have 
fixed before. If A is reduced the two intersection points 
diverge while for amplified A the intersection points ap¬ 
proach one another, i.e. their eccentricities change. This 
means that infinitely many physically distinct pairs of 








9 


orbits have the same frequencies as the originally fixed 
one if we allow the cosmological constant to take values 
in a certain interval. 

In the next section we discuss what happens to the 
degeneracy features when a spin of the particle is added. 

3.2. Spinning particles in the equatorial plane 


nominators have no zeros in the integration interval and 

V^{x) = (1 — + kix)(l + k2x){l + k3x){l + k^^x) 

(53) 

X (1 + fc 5 x)(l + k^x) (x — (a/ + ibi)) {x — {ai — ibj)) 

(54) 

with 


As before, we are interested in the influence of a parti¬ 
cle’s spin on its motion and dynamical properties. From 
(Hgi-® we obtain 


fj /”■“ 

y J Vp 

\/rPg (r) 

(46) 


\/rPg (r) 

(47) 



(48) 


with Pgir) from (|M 1 ) and 


Vtir) 




[g(r3 + (A,3_,),2)^(A,3_ 

(_A,3+,_2) (i+a, 2 ) 

^2 S^) 

(r3 + ( A^3 _ 1) (1 + Ag2) ■ 


(49) 

(50) 


The order of the polynomial rPg (r) under the root in the 
denominator of the integrand is 10. Therefore we now 
have a considerably more difficult kind of hyperelliptic 
integral than in the spinless case. 

Following the same procedure as in the non-spinning 

case, we transform the integrals with r = , , ^ - to 

obtain the resulting frequencies as functions of (p,e). 
Therewith, we are able to reparametrize the radial fre¬ 
quency, again, as a function of (fl^, e) so that the com¬ 
parison of different orbits and their frequencies is simpli¬ 
fied. Since the inclusion of the spin does not make the 
system easier to be solved, we encounter the same nu¬ 
merical problems as before, i.e. the numerical inversion 
to p(n 0 , e) and the divergencies close to the separatrices, 
as we have in the non-spinning case. Luckily, these prob¬ 
lems can be tackled also by the same method, only the 
expressions become more complicated. 

Employing Sochnev’s approach of approximating hy¬ 
perelliptic integrals we rewrite the integrals of (I46p and 
(l47ll with the substitution given in (|4T]) to obtain 


Tr.=A^rp,ra,A,s) f dx , (51) 

J — 1 yj Vj. 

Acj)=B^{rp,ra,A,s) / dx, (52) 


where V)® (x) and (x) are rational functions whose de- 


Tg - Tp 
Tg+Tp' 

_ Tg - Tp _ 

{rg - n) -I- (rp - ri) ’ 
rg - Tp 

{r2 - Tg) -I- {rg -P2) ’ 

_ rg - Tp _ 

{rg - ra) + {rp - Ta) ’ 

_ rg - Tp _ 

{rg - rp) + {rp - rp) ’ 

_ rg - rp _ 

{rg - r^) + {rp - r^) ’ 

which satisfy 0 < kg < k^ < k^ < ki < kg < ^ and 
— 1 < fcs < 0. As in the non-spinning case, we use 
Sochnev’s method for evaluating these integrals near the 
separatrices, see the Appendix. 

In Fig. [6] plots of the isofrequency pairing phenomenon 
for particles with two different spin values (0.1,—0.1) 
moving in the same Schwarzschild-de Sitter spacetime 
with A = 0.0005 are shown. The qualitative shape of 
the region of bound motion is similar to that of the non¬ 
spinning particle. The quantitative differences become 
apparent in the characteristic features of bound motion 
and its boundaries. Some basic facts that we have al¬ 
ready seen in the {p, e)-diagram are obvious, such as the 
shifts of the maximal eccentricity, of the ISCO and of 
the OSCO. However, what interests us is the impact on 
the domain where isofrequency pairing occurs and the 
question of whether further degeneracies due to the spin 
emerge. 

First, we qualitatively compare the size of the region 
where isofrequent orbits exist. It can be characterized by 
the azimuthal frequency at the ISCO (OSCO) and the 
one at the intersection of the COD line with the hori¬ 
zontal axis. Only orbits having an azimuthal frequency 
within this range do have isofrequent partners. The size 
of the allowed frequency interval decreases if the spin is 
chosen to be negative and increases if the spin value is 
positive. Although the region shrinks for negative spin it 
will never completely vanish as long as bound motion ex¬ 
ists. Therefore, the spin does not destroy this degeneracy 
in the fundamental frequencies of the orbital motion. 

Next, we compare the evolution of the maximal radial 
frequency for different spin values in Fig. [T] We notice 
immediately the shift of the entire curve closer to the 
center for positive spin and further away for negative 
spin values. This coincides with the trend of the shifts of 


fci = 

kg = 

ks = - 

ki = 
k5 = 
ke = 
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e 



e 



Figure 6: The figure depicts the phenomenon of isofrequency 
pairing for bound orbits in the e)-plane for a fixed 
A = 0.0005. The upper figure shows the chracteristics for 
s = 0.1 and the lower one corresponds to s = —0.1. The 
thick (red) boundary lines correspond to the separatrices and 
confine the region of bound motion. The (blue) solid lines 
inside these regions coorespond to constant values of Qr- The 
(green) dashed line represents the singular curve, i.e. the locus 
where the Jacobian determinant of the transformation from 
(p,e) to (fir, 11,^) vanishes. The Circular Orbit Duals (COD) 
marking the boundaries of the domains where isofrequency 
pairing occur are shown by the black dotted lines. 


the ISCO and OSCO. We also see the spin dependence of 
the value for for A = 0 (bold dots). More generally, 
if we choose a value for on the vertical axis and 

determine the intersection points with the three curves 
not only the position differs but also the corresponding 
A is not the same for the different systems. 

We will now investigate if the spin induces a further 
degeneracy in the sense that isofrequency pairs with dif¬ 
ferent spin values may have the same frequencies. We 
can deduce this from Fig. [6l Choose a line of constant 
radial frequency that lies in the isofrequency range. Then 
draw a line of constant azimuthal frequency that corre¬ 
sponds to a circular orbit located outside the range be¬ 
tween the ISCO and the OSCO and mark the two inter¬ 
section points. We know from the analysis of that 

the maximal radial frequency increases with positive spin 
and decreases with negative spin, keeping the cosmolog¬ 


^max 



Figure 7: The figure shows the change of 17™®’' in value as 
well as in position from the center when the value of A is 
varied. The (red) dots mark the Schwarzschild case with A = 
0. Along each curve A increases up to its maximal value when 
17™®’' vanishes. The black solid curve correponds to s = 0, the 
(green) dashed one to s = 0.1 and the (blue) dotted one to 
s = -0.1. 


ical constant fixed. Now, imagine the line of constant 
fir shrinking and spreading slowly for small negative and 
postive values for the spin parameter, respectively. The 
two intersection points will diverge for positive spin and 
approach one another for negative spin leading to six 
different orbits that have different eccentricities. For ex¬ 
ample, if the inital intersection points belong to a non¬ 
spinning particle and the spin is slightly changed to pos¬ 
itive or negative values, there exist infinitely many phys¬ 
ically distinct pairs of orbits having the same frequen¬ 
cies but different spin values. Even the spin direction 
is different in this example. If we allow both the cos¬ 
mological constant and the spin parameter to vary we 
get two-parameter families of isofrequency pairs with the 
same frequecies, parametrized by (A,s). 


4. CONCLUSIONS AND OUTLOOK 

In this work we investigated the characteristics of 
isofrequency pairing for both geodesic motion of a test- 
particle and the non-geodesic motion of a spinning test- 
particle moving in the equatorial plane of Schwarzschild- 
de Sitter spacetime. In contrast to the case without a 
cosmological constant, there exist two regions in the do¬ 
main of bound orbits where isofrequency pairing occurs. 
More precisely, it is not only the strong field regime that 
exhibits such a feature but also a region close to the out¬ 
ermost stable circular orbit. This is associated with a 
greater variety of zoom-whirl orbits than in the case with¬ 
out a cosmological constant: There are now not only or¬ 
bits that whirl near the periastron but also orbits that 
whirl near the apastron. 

Generally, adding a cosmological constant and/or the 
spin leads to the emergence of additional degeneracies in 
the frequencies. This occurs already for arbitrarily small 
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values of the cosmological constant and the spin. At least 
in principle, this additional degeneracy is of relevance in 
view of gravitational wave data analysis. Here we may 
think of an EMRI which, among all possible gravitational 
wave sources, is the closest physical realization of the 
dynamical system considered in this paper. Whenever 
isofrequency pairing occurs, knowledge of the fundamen¬ 
tal frequencies alone does not determine the shape of the 
orbit, i.e., additional information on the spectrum has to 
be taken into account. 

The most obvious plan for follow-up work would be 
to consider more general spacetimes, such as the Kerr- 
deSitter-NUT... spacetime. In addition, there are several 
other avenues for future studies on isofrequency pairing 
which we would like to mention briefly. 

From a theoretical point of view isofrequency pair¬ 
ing is of relevance in view of perturbation techniques. 
For example, in order to use KAM theory for perturbed 
integrable systems, certain non-degeneracy conditions 
have to be satisfied. The simplest version of these non¬ 
degeneracy conditions is obviously violated if there is a 
degeneracy in the frequencies such that other, more com¬ 
plex or more restrictive, conditions have to be tested. 
To mention another example, the feature of isofrequency 
pairing and the occurrence of a singular curve can be 
used to compare different approaches to the general rel¬ 
ativistic two-body problem, as it was already mentioned 
in 0. Methods such as the effective one-body approach 
or the post-Newtonian approximation can profit from the 
isofrequency pairing and its related characteristics. 

Also from a theoretical point of view, it is an interest¬ 
ing question to ask if there are spacetimes where three 
or more orbits with the same frequencies exist. In all 
examples treated so far there are only isofrequency pairs. 
(Here we are refering to the situation that all the pa¬ 
rameters of the dynamical system have been fixed which, 
for the cases treated in this paper, means fixing A and 
s.) As an attempt to find a candidate for isofrequency 
triples one could start with a Bertrand spacetime and 
perturb it a little bit. Bertrand spacetimes, which were 
introduced in [4l|, are spherically symmetric and static 
spacetimes in which the ratio of the radial frequency and 
the azimuthal frequency is a constant rational number q 
for all bound orbits, so they show the same total degen¬ 
eracy of the frequencies as the Kepler problem but now 
with q ^ 1. We are planning to search for isofrequency 
triples etc. in future work. 
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Appendix 

Approximation of hyperelliptic Integrals 


In this section we briefly describe how we use the ap¬ 
proximation scheme for hyperelliptic integrals developed 
by Sochnev [s^. He based the method on the approxi¬ 
mation of irrational functions by rational ones. In par¬ 
ticular, the irrational function c = ^^ciC 2 ...Cm can be 
approximated by the sequences {a„} and {bn} which are 
defined iteratively by 


ai 


Cl -f C2 -f ... -f Cm 
m 


and 


Cln+l 


(m - l)an -I- bn 

m 


bi 


CiC2‘“Cm 


a 


m— 

1 


1 ’ 


(55) 


^n+1 — 


-,m —Iz- 


a 


m — 1 
n+1 


(56) 


for n > 1. While {a„} approaches c from above, {&„} 
comes from below, i.e. oi > 02 >...>«„> c > > 

... > &2 > ^ 1 , and their common limit for n —> 00 is c. 

Let us consider a specific example related to our prob¬ 
lem of hyperelliptic integrals. The irrational function 
Vl + kx with |fc| < 1 is finite for —1 < x < 1. (These 
boundaries become important when we consider the 
boundaries of the integrals.) Therefore, we can approx¬ 
imate this function by choosing "^1 + kx = !^ciC 2 ...Cm 
with Cl = 1 + kx and C 2 = C 3 = ... = Cm = 1 and eval¬ 
uating the sequences {a„} and {bn} within the defined 
range of x. In this way we approximate our irrational 
function by rational ones. 

This method can be used for evaluating our hyperel¬ 
liptic integrals (l42)l . (|43ll where m = 2. To that end we 
have to approximate the function ^JVr{x) where Vr(x) is 
given by (1441) . If the absolute values of the k factors are 
far away from one, the procedure goes like this: First, we 
extract the factor (1 — x^) out of the radicand, arrange 
the k factors in decreasing order in their absolute values 
and group positive and negative fcs. Then we form sub¬ 
groups with m = 2 elements within each group where 
we supplement a factor of one when there are less than 
m = 2 elements in the subgroups, i.e. 


\/Vr{x) = \J\- x"^ (57) 

X y^(I — k 3 x){l + kix) v^(I -I- k 2 x) ■ I v^(I -I- k^x) ■ I 


Now we are able to compute the approximating sequences 
{an} and {bn} for each subgroup up to arbitrarily high 
order in n, always resulting in a rational function. Con¬ 
sequently, the integral that has to be solved can be ap¬ 
proximated by integrals of the form 


J R (^x, \/l — dx 


where R (x, Vl — x^) denotes a rational function of x and 
VT — x^. Using any of the three Euler substitutions or 
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elementary transformations which rearrange the form of 
the integral into tabulated ones, it is possible to solve the 
integral in terms of elementary functions. 

This gives a good approximation scheme as long as the 
absolute values of all ks are far away from one. However, 
we are interested in frequencies close to the separatrices 
corresponding to absolute values close to one for either 
or ki. Luckily, only a few modifications to the procedure 
are necessary to adapt it to this case. 

To begin with the integrals in (l42ll and (l43ll are divided 
into two integrals, where one runs from —1 to 0 and the 
other from 0 to 1. Then, it is not the factor (1 — that 
is extracted. Consider 


(1 — x)(l + x)(l + fc 3 x)(l + kix){l + k 2 x){l + k 4 x) 


where the k factors are already arranged in decreasing 
order in their absolute values and (1 — x^) is rewritten 
as (1 + x)(l — x). In the first integral (from —1 to 0) 
the product (1 +x)(l + fcsx) is taken out where fca is the 
greatest of the positive coefficients h leading to 



dx 

y'(l + x)(l + ksx) A{x, fci, ^ 2 , ki) 


where A{x,ki,k 2 ,k 3 ) = ■\/{l + kix)(l + k 2 x){l — a;)(l + k 4 x) 
has to be approximated by the same procedure as ex¬ 
plained above. The second integral (from 0 to 1) is rear¬ 
ranged in such a way that it yields 

1"^ dx 

Jo -^(1 — x)(l -I- k 4 x) S(x, fci, ^ 2 , fca) 

with B{x, ki,k 2 , ks) = {1 + x){l + fc 3 x)(l -I- fcix)(l -I- k 2 x). 

Here, the factor (1 -|- k 4 x) is pulled out together with 
(1 — x) because k 4 is the greatest of the negative coeffi¬ 
cients ki- Again, the remaining function is approximated 
resulting in a rational function in x. Therefore we obtain 
for our integral an approximation of the form 

Ri{x)dx ^ 1“^ i? 2 (x)dx 
i-i \/(l -I- x)(l -I- ksx) Jo yj{l- x)(l -I- k 4 x) 


the solution of which can be found in terms of elementary 
functions. In order to simplify the calculations further we 
may apply partial fraction decompositions to each of the 
rational functions providing us with integrals of the form 


dx 


-1 (I -I- a(fci, ^ 2 , k 4 )x)yj{l + x)(I -I- fcax) 
dx 

0 (I -I- 6(fci, ^ 2 , ks)x)^{l - x)(I -I- k 4 x) 


In the spinning case, we have to evaluate the integrals 
m and (15^ . i.e., we have to approximate the function 
y/V^{x) with I^®(x) given by (15^ . It is convenient to 
treat the two non-real zeros and the real zeros separately. 
Since 


\/(x - (a/ -I- ibj)) (x - (o/ - ibi)) 
is a real-valued irrational function it is possible to ap¬ 
proximate it by Sochnev’s method leading to 


\/(x - {ai -I- ihj)) (x - (a/ - ihi)) -A 



in first approximation, which proves to be of sufficient ac¬ 
curacy for our purposes. Therewith, the remaining terms 


— X ^){1 + kix){l + k2x){l + k3x){l + fe4x)(l -I- k5x)(l + kex) 

can be approximated analogously to the non-spinning 
case. Here, close to the separatrices either k 2 or —k^ 
approaches 1, and the approximation only contains ele¬ 
mentary integrals which, after partial fraction decompo¬ 
sition, reduce to integrals of the form 



dx 

(I -I- as(fci, ^ 2 , k4)x)sj{l + x)(I -I- k 2 x) ’ 
dx 

(1 -I- bs{ki,k 2 , k 3 )x)y/{l - x)(I -I- fcax) 


(58) 
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